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We study the effect of bias voltage on the nuclear spin polarization of a ballistic wire, which contains elec¬ 
trons and nuclei interacting via hyperfine interaction. In equilibrium, the localized nuclear spins are helically 
polarized due to the electron-mediated Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction. Focusing here on 
non-equilibrium, we find that an applied bias voltage induces a uniform polarization, from both helically po¬ 
larized and unpolarized spins available for spin flips. Once a macroscopic uniform polarization in the nuclei is 
established, the nuclear spin helix rotates with frequency proportional to the uniform polarization. The uniform 
nuclear spin polarization monotonically increases as a function of both voltage and temperature, reflecting a 
thermal activation behavior. Our predictions offer specific ways to test experimentally the presence of a nuclear 
spin helix polarization in semiconducting quantum wires. 

PACS numbers: 62.23.Hj, 75.75.-c, 73.21.-b, 31.30.Gs 


I. INTRODUCTION 

Magnetic structures are promising platforms for many mod¬ 
ern devices, e.g. memory,® sensors, 21 and quantum compu¬ 
tation hardware. 3 The opportunities to get an ordered mag¬ 
netic phase in the bulk and low-dimensional systems due to 
Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction 4 were 
studied in a number of theoretical and experimental works. 8 14 1 
The prominent feature of RKKY interaction in ID systems is 
the ordering of localized spins into a helix. 15 16 

When the current is driven through the system of electrons 
and nuclei, the spin polarization can be swapped between 
the two subsystems through the hyperfine interaction, lead¬ 
ing to dynamic nuclear polarization effects. 17 22 If the polar¬ 
ization of current carrying electrons and localized spins dif¬ 
fer, the spin-transfer torque arises pi 24 important for dynam¬ 
ics of domain walls 25 2 ^ and enhancing the tilting of the spiral 
structure in helimagnets. 27 Closely related is the dynamic nu¬ 
clear polarization, arising e.g. in helical edge states of topo¬ 
logical insulator. The backscattering of helical electrons can 
be of different origins, such as assisted by phonons, 28 mag¬ 
netic impurities, 29 or absence of axial spin symmetry. 30 It 
was shown that nuclear-assisted backscattering of electrons 
due to hyperfine interaction induces nuclear polarization when 
the current is driven through the edge states of topological 
insulator !^— I 

The main motivation for our work comes from the recent 
experiment by Scheller et al. , 33 where the conductance of a 
cleaved edge overgrowth GaAs quantum wire was measured. 
The measurements showed that the conductance of the first 
mode becomes e 2 /h at low temperatures instead of the naively 
expected 2 e 2 /h. This suggests the lifting of electron spin de¬ 
generacy. The possible explanation is the presence of a helical 
nuclear spin polarization that gaps out one subband and thus 
provides an electron spin selection. Further ways to confirm 
the presence of the nuclear spin helix were suggested theoreti¬ 
cally, by means of nuclear magnetic resonance, 34 nuclear spin 
relaxation, 35 and quantum Hall effect anisotropies 

In this work we propose and study a complementary 
method to detect nuclear spin helical polarization in the wire. 


It is based on the effect of bias voltage applied to the wire and 
therefore straightforward to perform experimentally. We in¬ 
vestigate how the bias voltage applied to the wire affects its 
nuclear spin polarization. We assume that at zero bias and 
finite temperature, nuclear spins are partially polarized into 
a helix due to the RKKY interaction. We find that an applied 
voltage induces a uniform nuclear polarization from both heli¬ 
cal and non-polarized nuclear spins available for nuclear spin 
flips via electrons. Therefore, upon increasing the voltage the 
helical nuclear polarization drops, while the uniform polar¬ 
ization grows, and the total polarization grows too. For small 
voltages and increasing temperature, the uniform polarization 
grows because of thermal activation of electrons, while the 
helical polarization dramatically drops in magnitude. Once a 
macroscopic uniform polarization has developed, the remain¬ 
ing nuclear spin helix rotates as a whole around the axis along 
the uniform polarization. Since the helical polarization affects 
the conductance of such systems p l 16 l 33 * 36 l these predicted fea¬ 
tures are expected to show up in the voltage and tempera¬ 
ture dependence of the transport currrent and thus they can 
be tested experimentally. Recently, cantilever-based magnetic 
sensing techniques have been reported which enable nuclear 
spin magnetometry of nanoscale objects such as the nano wires 
considered hereP 7 Such powerful techniques offer promising 
perspectives for direct experimental tests of the results ob¬ 
tained in this work. 

The paper is organized as follows. In Sec. [II] we present the 
Hamiltonian of our model. In Sec. Ill we describe the proper¬ 
ties of the electron bath. The derivation of the Bloch equation 
for the total nuclear spin in the wire is discussed in Sec. [IV] 
The resulting nuclear spin polarization and its dependence on 
the parameters of the system are presented and discussed in 
Sec. El Our conclusions follow in Sec. ED Additional infor¬ 
mation about our calculation is given in Appendix. 


II. THE MODEL 


We consider a one-dimensional electron gas and localized 
spins in a semiconductor nanowire. We will refer to these lo- 
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calized spins as nuclear spin in the following, however, they 
can be also of other origins, such as e.g. magnetic impuri¬ 
ties, etc. The electrons and nuclei interact via the hyperfine 
interaction described by the Hamiltonian 

H hyP = ^Apo 1 \^±(R±)\ 2 S(r - R)a ■ I , (1) 

where A is a hyperfine constant of the material, po is the 
nuclear spin density, 'ipx is the transverse part of electron 
wavefunction, r denotes the electron position along the wire, 
(i?, R±) is the position of the nucleus along the wire and in 
the transverse direction respectively, cr is an electron spin op¬ 
erator, and I is a nuclear spin operator (in units of h) with 
the magnitude /. We assume that the transverse part of the 
electron wavefunction 'ip±(R_\_) is constant in the wire cross- 
section, \'ip±(R_\_)\ 2 = 1/C, where C is the wire crossection 
area. We parametrize it alternatively by the number of nu¬ 
clear spins in the cross-section, N± = Capo, with a being 
the lattice constant. In GaAs po = 8/a 3 , a = 0.565 nm, 
A = 90 /ieV, / = 3/2, and N± is typically of the order of 
10 3 . Finally, we introduce N = L/a with L the wire length 
(typically of order microns), which gives NN± as the total 
number of nuclear spins in the wire. 

The total Hamiltonian reads 

H tot = -^d 2 r +H hyp , (2) 

where m is electron effective mass and Ti is the Planck con¬ 
stant. If the hyperfine interaction, Eq. (|T]), is weak on the 
energy scale of the electrons, its effects can be treated per- 
turbatively. The condition is quantified by A <C £p, where 
£p is the Fermi energy of the electron system. This condition 
is well satisfied in the cases we consider here. A Schrieffer- 
Wolff transformation on H tot perturbatively in H} iyp , i.e., in 
order A/ep, results to leading order in an effective interaction 
between the localized spins, the RKKY interaction!- 4 7 1138 

Hrkky = ^ ^ Ii ’ Jij Ij • ( 5 ) 

hJ 

Here, the indexes i. j label the nuclear spins and the RKKY 
coupling Jij = J( | Ri — Rj |) is related to the static spin sus¬ 
ceptibility of electrons (see Eq. (Cl) and below in Ref. 36), 
giving rise to the spatially dependent RKKY interaction. 

Let us rewrite Eq. 0 in the momentum represen¬ 
tation, defined through the Fourier transforms J q = 
exp [—iq(Ri — Rj)\Jij, with Ri E a, 2a,..., Na, and 
I q = exp(igi2i)Ii, with i E 1,..., NN±, and in both 
cases q E (27 t/N) x {0,1,..., N — 1}. We get 

Hrkky = -jy ^ ] Ig ' Jq I—g • (4) 

Q 

In one dimension, the RKKY coupling J q has a sharp mini¬ 
mum at momentum q = ±2kp, with kp = y/2m£p/h the 
electron Fermi wavevector. 15 16 Consider an approximation in 
which we neglect all values of J q with respect to the large 
(negative) value at this minimum, 

Hrkky — ~^J2k F (^2/c F ’ I— 2k F +' 1—2 k F ’ 1-2k F ) • (5) 



FIG. 1. A sketch of a conducting wire (yellow cylinder) with itiner¬ 
ant electrons (not shown) that couple to localized nuclear spins (red 
arrows) via hyperfine interaction. As a result, a helical nuclear po¬ 
larization emerges below a critical temperature. The blue spiral is a 
guide to the eye showing the direction of the helical polarization. The 
helical plane is chosen to be perpendicular to the wire axis (which 
need not be the case in general). 


To understand the spectrum of this Hamiltonian, we introduce 
linearly transformed spin operators, 

li = H u ^k F RiIii (6) 

with the matrix corresponding to a rotation by angle 0 
around a unit vector u. Inserting Eq. ([6]) into Eq. ([5} we get 

H rkky — J2k F (-^=0 * Jq =0 + 5 ( 7 ) 

where we define the vector components along u as I u = I • u, 
and perpendicular to it as I 1 - = I — I u u, and we separated 
the terms bilinear in the spin operators at finite momenta, 

B = Y [tift-q+Iiq-I-2 q + is Sn(q)(i2 q x£-2 g )- U }- 

q=±2k F 

( 8 ) 

The first term in the bracket of Eq. 0 describes the en¬ 
ergy of ferromagnetically coupled spins Ip. a configuration 
in which all these spins are collinear, along a vector per¬ 
pendicular to u, gives a minimal possible energy, of value 
NN]_J 2 k F I 2 - This configuration corresponds to a classical 
ground state of Eq. as well, as it saturates the energy lower 
bound obtained using the Parseval’s identity ^2 q \I q \ 2 = 
NN± JY | J/ 2 . Going back to the laboratory frame according 
to Eq. ([6]), the ground state corresponds to a helical configura¬ 
tion where the nuclear spins are oriented parallel to each other 
in the wire cross-section, along a direction which rotates in a 
fixed plane as one moves along the wire (for illustration, see 
Fig. [I}. We shall refer to this plane as the helical plane, with 
u being its normal unit vector. A unit vector h _L u gives 
the direction of the polarization within this plane at position 
R = 0. 

The finite momenta components, Eq. ([8]), arise from the 
choices of a definite helicity and the vector u in Eq. ([6]), which 
break the full spin rotational symmetry of Eq. 0. Namely, 
choosing a frame with helicity opposite to the ground state 
helicity would lead to a swap of the roles of i/p and J± 4fe . 
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Further, configurations where both helicities are populated 
lead to a lower energy gain. For example, choosing both 
with the same weight, gives in the laboratory frame a spin- 
density wave, i.e., a cos-like oscillation along a fixed vector, 
Ii = h cos (2 kpRi), which gives only half of the energy gain 
of a helical order. Such oscillating, rather than rotating, con¬ 
figuration corresponds to the first term in Eq. ([8]). To conclude, 
up to the spin rotational symmetry, which allows for arbitrary 
directions of u , and /- 1 , the ground state with ferromagneti- 
cally aligned Ii (helically ordered Ii) is unique. 

If the order is established, the expectation value of 0 is 
macroscopic, and we parametrize it by a polarization ph , 

<4=0> = NN ± Ip h h, (9) 

so that Ph = 1 corresponds to a completely ordered state. 
With this we reduce Eq. by the mean field approximation 
to a Hamiltonian describing a set of non-interacting spins 

h rkky ~j2^ b ?-I" ao) 

in the presence of the position-dependent internal field 

UnB™ = 2phN±I J2k F ^u,2k F Rih. ( 11 ) 

This concludes a simplified derivation of the reduction of the 
RKKY Hamiltonian, Eq. 0, into a set of non-interacting 
spins, Eq. ( p~Q] >, in an effective (mean) field, Eq. ( fTT) . A de¬ 
tailed analysis of the applicability of such an approximation 
was given in Ref. 36] based on the derivation of the spectrum 
of the full Hamiltonian Eq. 0, without employing a mean 
field ansatz. There it was found that this approximation, in 
essence neglecting the long wavelength magnons, is well jus¬ 
tified for sub-Kelvin temperatures and wire lengths relevant 
for mesoscopic experiments. 

As we consider the limit A <C £p, we adopt the Bohr- 
Oppenheimer approximation, assuming that electrons react 
instantaneously to the changes in nuclear spin subsystem. 
Consequently, we can consider the effect of the nuclear po¬ 
larization on electrons as an Overhauser fielcP 

Me B °v = 2 ^- £ 5(r - (12) 

3 

where fi e is an electron magnetic moment. Thus, the electron 
Hamiltonian is 

h 2 

H ei = + VeBov • er. (13) 

2m 

In Eq. ( [13] ) we do not include electron-electron interactions 
explicitly. In the following to evaluate the internal field 
we use Eqs. (C4), and (C5) from Ref. [36] In these equa¬ 
tions electron-electron interaction is significant (for example, 
for the critical temperature of the helical polarization 15 16 l 34 l 36 l) 
and therefore is included. 

To describe the nuclear polarization in the wire when a bias 
voltage is applied, we will first investigate the behaviour of 
one nuclear spin placed in an effective field of all others, Eq. 

and interacting with the bath of electrons described by 
Eq. ([13). 


III. HELICAL ELECTRONS AND FINITE VOLTAGE 


To find how a nuclear spin is affected by the electrons when 
the bias voltage is applied, we first consider the properties of 
the electron bath in the wire. As already mentioned in Sec. [IIJ 
the electrons are moving in the Overhauser field produced by 
the nuclear spins [see Eqs. <[13)]. As the nuclear spins 
form a helix in equilibrium, this particular Overhauser field, 
denoted by Bh , is also helical. Consequently, the electron 
spectrum is 


= h (k +k F ) ± p Jrr^ix\B\ + (14) 

2 m rri v 


where k is the electron wavevector, and £_ and £ + denote the 
lower and upper subbands respectively. They are split by the 
gap 2p e Bh at k = 0. The corresponding wavefunctions are 


= A=[e~ ik - r cos | |t) + e ikpr sin |t)], (15) 
«Vh(0 = A=[e ik * r cos |t) - e~ ik * r sin ^ |t)], (16) 


where cos Ok 

-mheBh 


h 2 kk F 


r-— -= and sin Ok = 

\J ( h 2 kk F ) 2 -\-(mfi e B h ) 2 

/ ,, 2 ,, and |t), |t) denote the spin states with 

v (h^kk F ) 2 +(mheB h ) 2 

spin up and spin down respectively, where u sets the quan¬ 
tization axis. These expressions of the wavefunctions can be 
simplified since typically the ratio A = PeBmax /£f 1, 
where B max is the maximum Overhauser field when all nu¬ 
clei are fully polarized along a given direction. For exam¬ 
ple, for a GaAs quantum wire /i e ^maai — 67 /xeV, while 
s F ~ 10 meV, which gives A 0.0067. Consequently, we 
can use A as a small parameter. 

We expand Eq. ( [15] ) in leading order of A and for the states 
within the partial gap we get 




J_r_ e i(fc-M r|£ fl ) ; A;>0; 

VL l e i ( k+kF ' )r la) ,k< 0, 


(17) 


where for right-moving electrons (k > 0) the spinor is |£r) = 
|t), and for left-moving (k < 0) it is |£l) = |t). Therefore, 
within our approximation the electronic states in the partial 
gap are helical: the spin is determined by the propagation 
direction, and is opposite for left-moving and right-moving 
electrons. 

Next, we consider the voltage applied to the wire and de¬ 
fine it as the difference between the chemical potentials for 
the left- and right-moving electrons (see Fig. [2]). Assuming a 
ballistic wire, the chemical potential of a given branch is con¬ 
stant in space. Concerning the nuclear spins, only the electron 
states from the partial gap are relevant, because we consider 
the case of small voltage, i.e. eV < 2p e Bh, where — e is 
the electron charge, and small temperatures T, i.e., k F T < 
2fi e Bh . This means that we adopt two approximations. First, 
we neglect the influence from electron states which are not 
in the partial gap (the upper (+) subband is neglected com¬ 
pletely), because their contribution to transport is exponen¬ 
tially small, proportional to exp[(— fi e Bh+eV/2)/ksT]. Sec¬ 
ond, we use Eq. ( [TT] ) for the electron wavefunctions, which 
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FIG. 2. Sketch of the energy spectrum given in Eq. GU and the 
direction of the electron spins in the presence of the helical Over- 
hauser field Bh and the uniform Overhauser field B u perpendicular 
to the plane of the helix. Red arrows denote the spin directions of the 
electrons in the lower subband, and the blue arrows label the spin di¬ 
rections for the upper subband. The coordinate system for the spins 
is formed by h and u shown in the right lower corner. The chemi¬ 
cal potentials for left- and right-movers are denoted as hl and hr, 
respectively. The voltage applied to the wire is gVrl = HR ~ Pl- 


means that we consider it in leading order of ksT/ep, eV/eF, 
and A. Therefore, for a description of the electron system 
in terms of a heat bath that causes the relaxation of the nu¬ 
clear spins, we take into account two branches: left- and right- 
moving electrons with spins \f) and |t), respectively. 

With the polarity as assumed in Fig. [2] the applied volt¬ 
age depletes the left (L) branch and increases the population 
of the right (R) branch. This imbalance in population opens 
up an additional phase space for the electrons to backscatter - 
predominantly from R to L. Because of the helical character 
of the states, such backscattering is accompanied by an elec¬ 
tron spin flip (from ||) to ||)). This, in turn, is enabled by 
the total spin-conserving hyperfine interaction Eq. {T}, so that 
each electron spin flip is compensated by a nuclear spin flip 
in the opposite direction. In this way a uniform nuclear polar¬ 
ization along the u direction is built up. We denote B u as the 
Overhauser field corresponding to this uniform polarization. 

The spectrum of electrons moving in the total Ovehauser 
field Bh + B u reads 




fi 2 (fc 2 + 4) ± 

2 m 



Ti 2 kkF 

m 



see Fig. [2] The asymmetry of the spectrum is due to the uni¬ 
form Overhauser field B u . The corrections to the wave func¬ 
tions, Eq. due to B u are negligible in leading order of 
A. 


We note that from Eqs. © and ( |T6| ) it follows that the 
electron spins become also polarized, thereby producing a 
Knight shift acting as an effective magnetic field B e - back on 


the nuclear spins. This Knight shift is defined as ( Hh yp ) e i = 
HnBj • Ij, where (...) e i denotes averaging over the eigen¬ 
states of the Hamiltonian in Eq. ( [13] ) with populations defined 
by the voltage. In this work, however, we can neglect B*j with 
respect to B^ produced by the RKKY interaction PEH 


IV. BLOCH EQUATION FOR THE TOTAL NUCLEAR 
SPIN IN THE WIRE 


To investigate the time-dynamics of the nuclear spins, we 
apply the standard Bloch-Redfield theory to our problem, 
which is valid for weak coupling between spin system and 
bath degrees of freedomp 9 ® as is the case here. First, we 
write down the Bloch equation for the average (Ij) of the j th 
nuclear spin. By applying Eqs. (7)-(ll) from Ref. |39]to our 
Eqs. and ( [13) , we get (for more details see Appendix 

0 


dt(Ij) = uj x {Ij) - Tjilj) + Tj, (19) 

where uj 3 = hnB^ /Ti determines the precession, the relax¬ 
ation tensor Tj the decay and the inhomogeneous vector term 
T j the stationary value of (Ij). Both, Tj and T j are ex¬ 
pressed in terms of time correlators (see App. A) 

i r°° 

JnM = —2 / e-^iSBn^B^^dt, ( 20 ) 

2n Jo 

where t is time, the indexes n, l label the components of 
the effective fluctuating internal field 5B defined via Hh yp — 
(Hhyp)el = SB ■ Ij. The time-dependence follows from the 
interaction representation 5B(t) = e ' lH eit/h^Be~ lHelt ) h . We 
note that above equations are valid for a spin 1/2. However, 
it is well-known 41 that the relaxation time of a spin into its 
stationary value does not depend on the spin length (in Born 
approximation). Thus, we will assume that our results apply 
for arbitrary spins. 

As follows from Secs. [II] an d|m| we can define the expec¬ 
tation value of a nuclear spin at position R = 0 as 

(I 0 ) /1 = p h h + p u u, (21) 

where 0 < Ph,u < 1 denote the polarizations along the two 
orthogonal directions h and u , respectively. 

We also introduce position-independent tensors r 0 and T 0 
in the rotated frame defined by the rotation matrix R) u 2 k F R ., 
via 


Tj — ^ U ,2k F R j ToR^ 2kF R., 
Tj = T^u,2k F RjTo. 


( 22 ) 

(23) 


Having Eqs. ©, ©, ©, and ( [23) , we can describe the 
time-evolution of the nuclear spin I j in the rotated frame. 

Eventually we are interested in the dynamics of the total 
(macroscopic) polarizations, rather than the one of an indi¬ 
vidual nuclear spin. We therefore introduce the total nuclear 
spin in the rotated frame V R) u 2kpR . (Ij) = NN_\_ (I 0 ), and 
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write the equation of motion for it using Eqs. <[TT]i, ([19]), 
and d23l. We get 


dt(Io) = -fi((/ 0 ) • h)h xu — r 0 (J 0 ) + T 0 , (24) 


where we denoted ft = Ip u \J2k F \/h • The first term implies a 
rotation of the helical direction h , around the axis n with fre¬ 
quency £1 This can be seen by introducing a time-dependent 
vector h(t) = 7where a(t) = Qdr. In the Born- 
Oppenheimer approximation, the tensors Tj and T j are func¬ 
tions of the instantaneous values of h and u, so we write 


r j = n uMt) r ( 26 ) 


With this the time evolution of (Jo) in the rotating frame, 
(J 0 ) = (Jo), is described by 


a t (Jo) = -f 0 (Jo) + T 0 . 


(27) 


Using this equation and properties of the electron bath dis¬ 


cussed in Sec. Ill we can describe the polarization of the 
nuclear spins in the wire as function of temperature and volt¬ 
age. 


J 


V. RESULTING POLARIZATIONS 


To find the polarizations ph and p u from Eq. ( [27] ) we now 
evaluate the tensors T o and To explicitly. For that we first 
evaluate the correlator J n i(uj). Using Eqs. (17) and (18) we 
get 




Ml 


ab 


nl 


32h 3 7rv 2 F Nl ^ Mnl Qab ’ 

b - L a,be{L,R} 

(28) 

(falOnl&X&Mfo), 

(29) 


Qab = J dsf (e + eVba/2) [1 -f(s + hco- eV ba /2 )] ,(30) 


where eVb a = Pb ~ P>a is th e difference between chemical 
potentials of branch b and a, with a and b denoting L (left- 
movers) or R (right-movers). Here we also use the Fermi 
distribution function f{e) = [exp [e/^ksT)] + l] -1 . As was 
mentioned in Sec. [In] we consider voltages and tempera¬ 
tures smaller than the partial gap 2 p e Bh given by the helical 
polarization. Therefore the term f{e + eV& a /2)[l — f{e + 
Tiuo — eVb a /2)] allows us to consider only the energy window 
of ±p e Bh around £p, because f(e) decays exponentially for 
e/ksT 1. Consequently, we approximate the electron den¬ 
sity of states (per spin) by v(e) ~ v(ep). Up to first order in 
A, we have v{ep) = 1/(tt hvp), where vp = ep/(hkp) is 
the Fermi velocity of the electrons. 

Having obtained it is straightforward to calculate 

r 0 and T 0 , using Eqs. ( |A2| )-( |A4| ) and ([28]) -([30]). We can then 
solve Eq. © for the steady state polarizations (keeping ujq 
as a constant) and obtain 


AhuJo 

(Hu o - eV) coth (^f?) + (hu>o + eV) coth + 2 hw 0 coth ’ 

_ (B<»o-eV)coth( ' i 3 °-° v )-(nh;o + ey) coth) _ y 

^ (hwo—eV) coth ^ t ^gT V ) +(fi<*>o+eV) coth ^ ) +2huo coth ^ ) 

(huJo - eV ) coth ( AfcsT 7 ) + (^o + eV) coth ) + 2eV coth 


(31) 


(32) 


However, from Eq. it follows that Tiujq = 

2 PhN±IJ 2 k F h, i.e., ujo depends on p This leads to non¬ 
linear algebraic equations for two unknowns, and ph , 
which we solve numerically using material parameters for 
GaAs (analytical expressions for small deviations of the po¬ 
larizations are given below). We plot the values obtained in 
this way and discuss their behaviour as a function of voltage 
and temperature, the experimental parameters that are most 
directly accessible. 

The voltage dependence of the polarizations is shown in 
Fig. [3] We can see that the polarization p u grows faster with 
voltage than pj h decay s, therefore the overall polarization of 
the nuclei \/p\ + p\ grows with voltage, too. This means 
that the nuclear spins are more polarized when a voltage is 
applied than when they are in equilibrium at the same tem¬ 


perature. We also note that having a non-zero component p u 
means that nuclear spins have a conical polarization, rather 
than a helical one. To plot Fig. [3] we used Eqs. (C4), and 
(C5) from Ref. [36] as was mentioned above, where the de¬ 
pendence of huo o on temperature is described in detail. To 
evaluate huo we used the characteristic values for GaAs: the 
Fermi velocity vp = 2.3 x 10 5 m/s, and the number of nuclei 
in the wire cross-section N± = 1300. For the expression for 
J 2 k F taken from Ref. [36] we use the electron-electron inter¬ 
action Luttinger liquid parameter K p = 0.2 and the absolute 
value of spin / = 3/2. For the constants described above and 
at T = 90 mK and p u = 0.1 the rotation frequency of the 
nuclear spin helix isU«1.5xl0 6 s -1 . 

It is natural to e xpect that high temperature destroys the nu¬ 
clear helical order. 15 ^ 34 36 Indeed, Fig. [ 4 ] shows that the heli- 
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FIG. 3. (a) The voltage dependence of the polarization ph along he¬ 
lical direction h (blue), polarization p u in the direction of u perpen- 
di cular to th e helix plane (red), and the overall polarization of nuclei 
\/Pu + Ph (black), (b) Enlarged from (a) the voltage dependence of 
Ph and y/pu + Ph- We use T — 90 mK and other parameters as 
given in the text. 
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FIG. 4. Plot of the temperature dependence of the polarization p u 
(upper panel, red) and the polarization ph (lower panel, blue). For 
these plots the same parameters were used as in Fig. [3]and the applied 
voltage is eV = 0.5 /xeV. We note that our calculation is valid for 
eV, ksT < 2 n e Bh, therefore the smallest value of ph we consider 
here is ph — 0.2. 


cal polarization ph decays with temperature and then drops in 
magnitude around T ^ 109 mK. As our calculation is valid 
for eV, ksT < 2 p e Bh, the smallest value of ph allowed by 
self-consistency for our parameters is ph — 0.2. From Fig. 
[4] it also follows that the polarization p u grows with temper¬ 
ature. This growth is explained by the fact that due to higher 
temperature the electron states with higher energy become oc¬ 
cupied. This makes the nuclear spin flip more probable. It is 
obvious that there is a temperature where the polarization p u 
gets destroyed, however, for the range of temperatures given 
in Fig. [4 \p u does grow, whereas the helical polarization ph de¬ 
cays significantly. The decay of ph with temperature is rapid, 
while the growth of p u is less pronounced. Therefore, the 
overall nuclear polarization in the wire strongly decays with 
increasing temperature. For the parameters we used for Fig. 
[4] the effect of temperature on ph is stronger than the one of 
a finite voltage. The initial temperature scaling of ph away 
from unity (see Fig. [4]) can be obtained readily from Eq. ( [3l] ) 
by treating 1 — ph as a small perturbation. This yields 


Ph 


2 

1 -^ 

1 e T9 

1 - 2e - ^, 


1 2T9 Se °h ( 2 T9 ) 


(33) 

(34) 


where the first equality holds well for the temperature inter¬ 
val 60mK < T < 90mK, while the second one is a good 
approximation for 60mK < T < 80mK. Here we denoted 
g = 3- - f AKp , and the temperature-independent param¬ 


eter k = 2N±I\ J<ik F \T 9 ~ 1 /kB depends on the material and 
geometrical properties of the sample (see Eq. (11) and Eqs. 
(C4) and (C5) of Ref. 36). For K p = 0.2 (chosen for the 
plots) we get g = 2.4 (we recall that K p = g = 1 corresponds 
to vanishing electron-electron interactions). 

The initial decrease of ph due to voltage in Fig. [3] for V < 
3 peV scales as 


Ph ^a-jV 2 , 


(35) 


where a and 7 depend on material and geometrical parameters 
of the nano wire and on temperature. 

Finally, we mention that recent progress in nuclear spin 
magnetometry on nanowires 37 has opened the perspective to 
measure the nuclear spin polarizations directly and thus to 
test the predictions made here. Moreover, due to the heli¬ 
cal nuclear polarization which acts on electrons as an Over- 
hauser field Bh there is a partial gap in the electron spectrum 
[see Eq. ©]. As a result, the conductance of a ballistic 
nanowire is less than 2 e 2 /h for sufficiently low temperatures 
and V < 2 /i e B h jHEMIIlHl As was shown above, the polariza¬ 
tion ph, and consequently Bh, decrease with increasing volt¬ 
age and temperature. We thus expect qualitatively that the 
conductance of the wire will increase with the decrease of the 
partial gap 2fi e B h oc p h W 
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VI. CONCLUSIONS 

We have shown that due to the hyperfine interaction be¬ 
tween electrons and nuclei in the wire the applied voltage 
changes the form of the nuclear polarization and its ampli¬ 
tude. Assuming that in equilibrium there is a helical nuclear 
polarization ph present in the wire due to RKKY interaction, a 
bias voltage induces a uniform polarization p u perpendicular 
to the helix plane. Due to this polarization the nuclear spin he¬ 
lix starts to rotate around the axis perpendicular to the helical 
plane. When a non-zero polarization p u buidls up, the nuclear 
polarization changes from helical to conical. 

We have also presented the voltage dependence of p u and 
Ph and seen that p u increases with voltage, whereas ph de¬ 
creases. Following from these two effects the overall nuclear 
polarization in the wire grows with voltage. Remarkably, p u 
grows with temperature in the considered range of temper¬ 
atures. This is because the nuclear spin flip becomes more 
probable as electrons occupy higher energy states. This ther¬ 
mal activation effect is noticeable for the considere d regime 
hujQ > eV. The growth of the overall polarization \Jp\ + p\ 
with voltage and the growth of p u with temperature are in¬ 
triguing and a priori non-obvious effects. The polarization 
effects predicted here might be observed in transport experi¬ 
ments 33 or more directly via cantilever based nanoscale mag- 
netometry. 37 
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Appendix A: Bloch equation for one nuclear spin 

To write down the Bloch equation for the total nuclear spin 
in the wire, we use Eqs. (7)-(ll) from Ref. [39j Here we 
present them adopted to our case of a nuclear spin interacting 
with the bath of electrons and placed into the effective field 


produced by all other nuclear spins in the wire. 

The Bloch equation for the nth nuclear spin reads 

9t{In) = x {In) ~ r n {I n ) + Y n . (Al) 

To express tensors T n and Y n we introduce a unit vector l 
along <jj n , i.e. u? n = uj n l. The tensor T n consists of a dephas¬ 
ing part Y d n which comes from energy conserving processes 
and a pure relaxation part T r n , which comes from the energy 
exchange with the bath 39 43 (played here by the electron sys¬ 
tem), 


r iij= o) - ( A2 ) 

1'n,ij = ~ ^p^q)^pq{^ri) — (A3) 

{$ip hlp)*?Tpj{^n) ^ij^kpq^k^pq^n) T~ 

^ipq Ip^Lqj (cc n ) ]. 

Here, the indexes z, j denote components of tensors, and we 
use the Einstein convention of summation over repeated in¬ 
dexes. Further, e pq k is the Levi-Civita symbol and Sij the Kro- 
necker delta, while Ik denotes the kth component of vector L 
The inhomogeneous part of the Bloch equation Y n read J 39 * 43 1 

n,i ~ 2 (^ n ) ^ipql'pqi^ n) T" (A4) 

+ e iqklklp[I'pq{bJn) ~ Ipq{®)]) ■> 

where i denotes the component of Y n . The terms 
Xf^iuS) are defined as 

M = R 4Jij M ± Jij (-w)], (A5) 

I^(w) = Im± Jij(-w)]. (A6) 

The term J, : j (u>) is the Laplace transformation of the correla¬ 
tor of the fluctuating fields 5B at different times, 


Ji 


i r 


= ^72 / (A7) 


where 5B(t ) = e lH eit/h^ e -xH el t/n^ u s j n g <A1) we 
expressed the Bloch equation for the total nuclear spin in the 
wire resulting in Eq. ([24|). 
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